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ABSTRACT 

Thermal  or  equilibrium  processes  contribute  to  the 
total  ionization  of  the  exhaust  plasmas  attending  a  missile 
in  flight.  Such  equilibrium  ionization  is  the  principal  con¬ 
tributor  to  the  total  ionization  in  some  cases  of  interest. 
Calculations  of  the  equilibrium  values  of  ionization  for  vari¬ 
ous  rocket-propellant  components  are  presented  in  graph¬ 
ical  form  to  make  numeric  values  readily  available  and  to 
simplify  their  utilization. 
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THERMAL  IONIZATION  OF  ROCKET  EXHAUST  PLASMAS 


The  plasmas  attending  a  missile  in  flight  may  stem  from  a  number  of  sources,  such 
as  from  shock  ionization,  chemi-ionization,  thermal  ionization,  photoionization,  etc.  One 
of  the  important  ionization  mechanisms  is  the  thermal  (or  equilibrium)  process.  Any  of 
the  ionization  sources  may  control  in  a  particular  case,  but  only  the  thermal  ionization 
will  be  presented  in  this  report. 


The  subject  of  equilibrium  excitation  and  ionization  is  considered  in  treatises  on 
physical  chemistry  (see  item  1  in  bibliography,  for  example).  It  has  been  treated  in  some 
detail  by  Saha  (2)  and,  for  this  reason,  is  sometimes  called  Saha  ionization.  For  the 
simple  case  of  a  single  ionizing  component,  the  expression  for  the  ionization  is  written  as 
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loRio  (Tl  -  6-49  +  lo^io 


where 


X  =  fraction  of  the  ionizable  component  which  is  ionized 
P  =  partial  pressure  of  the  ionizable  component  (atmospheres) 

I  =  ionization  potential  of  the  ionizable  component  (electron  volts) 

T=  temperature  (degrees  Kelvin) 
u,.  =  weight  factor  for  the  electron  (=2) 

=  weight  factor  for  the  ion 
a.,  =  weight  factor  for  the  atom 

and  the  weight  factor  is  2j  ^  1,  where  j  is  the  inner  quantum  of  the  ground  term  for  the 
particular  atom  or  ion. 


If  ttiero  is  more  than  one  component  that  is  ionized,  either  in  different  excited  states 
of  a  single  typo  of  atom  or  in  different  atoms,  the  relation  then  becomes 
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^  l<>R,o  (T)  -  6.49 


1  nuid 


where 

=  fraction  of  the  particular  ionizable  component  that  is  ionized 

i\  =  partial  pressure  of  the  particular  component 

f  =  fraction  of  free  electrons  to  the  original  number  of  undecomposed 
atoms. 

Similar  equations  are  required  for  each  excited  stati'  of  each  type  of  atom  nr  radical  con¬ 
tributing  to  the  total  ionization.  A  simultaneous  solution  of  the  set  of  equations  is  then 
required  for  a  solution.  Further  refinements  are  required  tor  multiple  ionization. 

A  knowledge  of  the  conditions  that  exist  is  usually  inadequate  to  permit  a  high  degree 
of  accuracy  in  the  calculation  of  chemical  kinetic  processes  wliich  lead  to  the  observed 
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ionization.  As  a  consequence,  the  first  of  the  above  equations  with  the  assumption  of  a 
single  source  of  electrons  is  adequate  for  most  computations.  Usually  the  accuracy 
desired  from  the  calculations  is  such  that  further  simplications  can  be  made.  As  a  con¬ 
sequence,  plotted  values  of  ionization,  such  as  those  illustrated  in  the  following  figures, 
are  frequently  of  utility. 

The  solution  of  the  equilibrium  equation  gives  curves  as  plotted  in  Fig.  1,"  which 
illustrates  the  steep  rise  in  ionization  at  the  lower  temperatures.  The  ionization  of  the 
total  gas  rather  than  of  the  ionizable  component  is  usually  of  more  interest.  The  curves 
of  Fig.  1  appear  as  shown  in  Fig.  2  when  expressed  in  terms  of  the  total  gas.  Use  of  an 
inverse  temperature  scale  simplifies  the  graphical  presentation,  as  shown  in  Figs.  3,  4, 

5,  and  6.  The  weight  function  lo6io(g^g;  /e.J,  which  has  a  value  of  zero  for  the  alkaline 
earth  metals  and  for  aluminum,  i.e.,  for  most  of  the  curves  presented  in  this  report,  was 
assumed  to  be  zero  for  the  general  curves.  Thus,  in  Figs.  5  and  6  and  other  similar  plots, 
a  nonzero  value  of  this  function  would  cause  a  relatively  small  shift  of  all  the  curves 
along  the  temperature  axis.  The  effect  is  illustrated  in  Fig.  7.  Similarly,  a  change  in 
the  pressure  causes  a  shift  in  the  T  direction  as  illustrated  in  Fig.  8. 

The  curves  for  the  ionization  of  all  gases  are  very  similar  when  plotted  as  a  function 
of  the  ionization  potential  divided  by  the  temperature  (l/  T).  This  is  illustrated  in  Fig.  9 
for  various  single -component  gases  having  widely  varying  ionization  potentials  (11.0  to 
2.0  ev)  and  for  a  total  pressure  of  one  atmosphere.  Usually,  the  range  of  ionization 
potentials  of  the  elements  of  interest  are  much  narrower  than  the  total  spread  illustrated 
in  this  figure  and  fall  within  the  range  of  the  two  inner  curves,  i.e.,  between  cesium 
(3.893  ev)  and  aluminum  (5.984  ev).  Thus,  a  single  intermediate  curve  would  give  results 
accurate  to  within  30  percent  for  most  cases.  (This  accuracy  is  sufficient  for  most  cases 
of  interest.)  A  set  of  curves  for  an  average  ionization  potential  is  illustrated  in  Fig.  10. 

K  results  at  pressures  other  than  one  atmosphere  are  desired,  the  curve  should  be  shifted 
horizontally  by  an  amount  determinable  from  Fig.  8. 

In  many  cases  of  interest,  the  electron  density  rather  than  the  fraction  of  electrons 
present  is  of  interest.  Curves  for  the  electron  density  at  twelve  atmospheres,  the  approxi¬ 
mate  pressure  at  the  throat  of  many  rocket  motors,  are  shown  in  Fig.  11  for  various  ion¬ 
izable  components  of  interest.  The  decrease  in  electron  density  at  high  temperatures  is 
a  consequence  of  the  reduced  gas  density  with  increasing  temperature  and  constant  pressure. 

The  ionization  in  a  gas  at  one  atmosphere  pressure  containing  cesium  is  shown  in 
Fig.  12.  The  curves  for  the  various  ionizing  components  may  be  combined  on  a  single 
family  of  curves  when  the  electron  density  is  plotted  against  the  dimensionalized  parameter 
I  T,  similar  to  those  used  for  the  fraction  ionized.  Curves  of  tliis  general  type  for  two 
pressures,  one  and  twelve  atmospheres,  are  shown  in  Fig.  13.  A  simpler  curve  with  suffi¬ 
cient  accuracy  fur  most  applications  is  illustrated  in  Fig.  14. 

Caution  must  be  exercised  in  applying  these  curves,  particularly  because  of  the  large 
number  of  cases  in  which  they  are  inapplicable  because  equilibrium  processes  are  not 
controlling.  Further  discussions  of  their  applicability  may  be  found  in  papers  referenced 
in  the  bibliography. 


'■.•Vll  figvirrs  .irc  located  at  the  end  of  this  report. 
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APPENDIX 


CALCULATION  OF  SINGLE  lONIZABLE  COMPONENT  X 


The  equilibrium  equation  for  a  single  ionizable  component  may  be  rewritten  in  the 
form 

j0-(6.49  +  5041  I/T) 
pz-j-s  2  +  io-(6.49  +  5041  I/T) 

For  p/t®  ^  £qj.  either  small  pressure  (p)  or  large  temperature 

(T),  this  equation  gives  x  -  1,  as  shown  on  the  various  curves  of  Figs.  1  to  10. 

Similarly,  for  T  small  (or  I  large)  so  that  10 «  p/T®  ^  and  5041 
I  /t  »  6.49, 

X  -  '*  10“(5f41  I  2 

Since  T  varies  slowly  compared  to  10'(®°'*'  ^  x  x  10"(2.S20  i  t)^  ^j.  ^  ^  -(2520  i/t). 

Thus,  on  the  log  linear  plots  used,  the  ionization  varies  approximately  linearly  with  l/T 
throughout  much  of  the  region,  as  shown  on  the  curves  using  an  inverse  temperature  scale. 
A  single  set  of  curves  will  suffice  for  all  ionization  potentials  whenever  this  approxima¬ 
tion  applies.  Figures  9  and  13  show  the  range  of  values  encountered  for  various  ioniza¬ 
tion  potentials.  Figures  10  and  14  show  appro.ximate  values  suitable  for  many  purposes. 

An  approximate  equation  may  also  be  written  in  the  form 

j0-(fi.49  *  5041  !.T)  '  ^ 

"  1_  . 

whenever  p  'T®  ^  '  ’T'.  Thus, 

xp'  ^  jQ-l  2(6.40  5041  I  T) 

T5  4  - 

or 


>orio  ’‘  =  -12  loRjn  P  +  5  4  loRp,  T  -  1  2(6.4')  ^  .S041  1  T). 

Thus,  a  change  in  pressure  p  shifts  the  value  of  the  ionization  much  in  the  same  manner 
as  a  negative  change  in  temperature.  Such  a  shift  is  illustrated  in  Fig.  8. 
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Fig.  1  -  P'rartion  of  ionizablc  roniponcnl  which  is  ionized 
varies  with  partial  pressure 

ri’iL'  f  rac  t  •.ii:;  ol  .i  Io'a  -  .  oi;;  z.i  t :  (jt.  -  jrot  t-n*  .a  i  l  o .  ].‘on  c  nl  'Ahieh  is  loiii/.i’fl 
.\t  a  t  rj-:'.  pt- rat  u  r  f  <it’C  r  <‘<i  s«- s  ,lv  th**  j)art;.il  prrssari'  of  thf'OiTi- 

;iu  Th<-  I’r.utior.  c>:  'Ahich  is  •.on;/.«  ri  as  a  {'ip.l- 

tion  o!  .*.  >  partial  pr  t* sur  «•  ara:  Mi*  *  ral'i  rr  ilhislraW’Ti  m  thf 

Miv.iiv  i)t  L-irvrs  sho’.M.  .^bovt-.  Krt-r.  tliomh  thr  traction  ;  fMi./’-d 
(ii’v' r  I’a  sr  <  'Ailh  i  n  c  a  s  .  n  c  pr  s •  i  r  r ,  tin-  total  '.‘iLLlron  (irnsity 
nc  r  ca  <  s  bfcaiM-  o:  th«*  in*.  asini;  cimsitv  of  parent  .tton.s. 


MOLE  FRACTION  ELECTRONS  IN  GAS 
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TEMPERATURE  (DEGREES  KELVIN) 

Fig.  2  -  Mole  fraction  of  total  gas  which  is  ionized 
increases  with  partial  pressure 

Thf  ciirvf.s  of  Fig.  1  arr  rcplottod  tin-  fraction  of  total  jjas  which 

IS  •.uni/’.cd  .IS  tin  ordinate.  Plots  such  as  this  v.  .Ill  the  elitiron  dcn- 
sitv  expressed  as  a  fraction,  or  sometimes  called  a  mole  fraction, 
are  of  particular  utility  when  considering  frozen  expansion  of  a  gas. 
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Fig.  3  -  Potassium  ionization  plotted  on  an  inverse  temperature  scale 


Fig.  5  -  Low  ionization  potentials  produce  plasmas  at  low  temperatures 

No  normal  t-lemi-nt  has  ionization  potentials  this  low.  However,  the  electron  affinity  of  some 
atoms  and  ions  produces  negative  ions  by  electron  capture,  such  ions  having  low  ionization 
potentials.  Theycapture  electrons  at  low  temperatures  and  give  them  up  at  highe  r  tempe  rature  s . 
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Fig.  10  -  A  family  of  partial  pressure  curves  suffices  for  most  calculations 

rhf  value  of  the  paranieters  contributing  to  thermal  ionization  are  not  no  rmal  ly  known  with  suf¬ 
ficient  accuracy  to  permit  accurate  calcuJation.s  of  the  electron  density.  In  such  cases,  a  single 
f.aniily  of  partial  pressure  curves  plotted  against  i/t,  using  a  mean  of  the  family  shown  in  Fig. 
8,  gives  result.s  with  an  accuracy  sufficient  for  most  purposes. 
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TEMPERATURE {THOUSANDS  OF  DEGREES  KELVIN) 


POTASSIUM  FRACTIONAL  TOTAL  PRESSURE-12  ATMOSPHERES 
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TEMPERATURE (THOUSANDS  OF  DEGREES  KELVIN) 

Fig.  11(c)  -  Electron  density  for  sodium  ionization 


TEMPERATURE (THOUSANDS  OF  DEGREES  KELVIN) 

Fig.  IKd)  -  Electron  density  for  aluminum  ionization 
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IONIZATION  POTENTIAL/ TEMPERATURE  (ELECTRON  VOLTS  PER  DEGREE  KELVIN) 

Fig.  13(b,i  -  Plasma  strength  for  various  ionization  potentials 
at  twelve  atmospheres  pressure 


